Ecology and Evolution 

Meta-analysis reveals complex marine biological responses 
to the interactive effects of ocean acidification and 
warming 

Ben P. Harvey 1 , Dylan Gwynn-Jones 1 & Pippa J. Moore 1,2 

1 1nstitute of Biological, Environmental, and Rural Sciences, Aberystwyth University, Aberystwyth, UK, SY23 3DA 
2 Centre for Marine Ecosystems Research, Edith Cowan University, Perth, Australia, 6027 



Keywords 

Climate change, meta-analysis, multiple 
stressors, ocean acidification, ocean 
warming, synergistic interactions. 

Correspondence 

Pippa J. Moore, Institute of Biological, 
Environmental and Rural Sciences, Edward 
Llwyd Building, Aberystwyth University, 
Aberystwyth, UK, SY23 3DA. Tel: +44 (0) 
1970 622293; Fax: +44 (0)1970 622350; 
E-mail: pim2@aber.ac.uk 

Funding Information 

Funding provided through a Marie Curie 
Career Integration Grant PCIG10-GA-201 1- 
303685 and an Institute of Biological, 
Environmental and Rural Sciences Ph.D. 
studentship. 

Received: 5 November 2012; Revised: 23 
January 2013; Accepted: 29 January 2013 

Ecology and Evolution 2013; 3(4): 1016-1030 



Abstract 

Ocean acidification and warming are considered two of the greatest threats to 
marine biodiversity, yet the combined effect of these stressors on marine organ- 
isms remains largely unclear. Using a meta-analytical approach, we assessed the 
biological responses of marine organisms to the effects of ocean acidification 
and warming in isolation and combination. As expected biological responses 
varied across taxonomic groups, life-history stages, and trophic levels, but 
importantly, combining stressors generally exhibited a stronger biological 
(either positive or negative) effect. Using a subset of orthogonal studies, we 
show that four of five of the biological responses measured (calcification, pho- 
tosynthesis, reproduction, and survival, but not growth) interacted synergisti- 
cally when warming and acidification were combined. The observed synergisms 
between interacting stressors suggest that care must be made in making infer- 
ences from single-stressor studies. Our findings clearly have implications for the 
development of adaptive management strategies particularly given that the fre- 
quency of stressors interacting in marine systems will be likely to intensify in 
the future. There is now an urgent need to move toward more robust, holistic, 
and ecologically realistic climate change experiments that incorporate interac- 
tions. Without them accurate predictions about the likely deleterious impacts 
to marine biodiversity and ecosystem functioning over the next century will not 
be possible. 
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Introduction 

The concentration of atmospheric carbon dioxide (C0 2 ) 
has increased from 280 ppm in preindustrial times to a 
present day level of 391 ppm (Le Que're' et al. 2012). Over 
the last 100 years, this has led to changes in global sea 
surface temperatures (+0.74°C) and ocean carbonate 
chemistry (Orr et al. 2005), which have included ocean 
acidification by 0.1 pH units (Caldeira and Wickett 2003; 
Kleypas et al. 2006). By the year 2100, sea surface temper- 
atures are expected to rise by a further 1— 4°C while 
increased C0 2 (aq) will result in the decreased availability 
of carbonate ions and a further reduction in pH by 0.3- 
0.5 units (Caldeira and Wickett 2005; IPCC 2007; Goo- 
ding et al. 2009). These changes in temperature and ocean 
carbonate chemistry are considered two of the greatest 



threats to marine biodiversity (Kleypas et al. 1999; Doney 
et al. 2009), leading to changes in the physiological per- 
formance of individual organisms, which will in turn alter 
biotic interactions, community structure, and ecosystem 
functioning. 

A range of marine biological responses have already 
been observed in response to ocean warming including 
hypoxia (Portner and Knust 2007), coral bleaching (Hoe- 
gh-Guldberg et al. 2007), species range shifts (Parmesan 
and Yohe 2003; Root et al. 2003), changes to phenology 
(Walther et al. 2002), and reduced organism body size 
(Daufresne et al. 2009). Experimental manipulations sim- 
ulating predicted future ocean temperatures have sug- 
gested that warming will also lead to increased metabolic 
costs for plants and animals (O'Connor et al. 2009), 
increased consumption rates (Sanford 1999), and changed 
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food-web structures (Petchey et al. 1999). Observed 
responses of marine organisms to recent ocean acidificat- 
ion are limited (but see Iglesias-Rodriguez et al. 2008b; 
Moy et al. 2009), but are expected to become increasingly 
apparent in the next 50-100 years (Doney et al. 2009; Fe- 
ely et al. 2009). Experimental evidence, however, suggests 
that responses are likely to be varied (Hendriks et al. 
2010; Kroeker et al. 2010) and will include hypercapnic 
suppression of metabolism (Christensen et al. 2011), acid 
-base balance disturbances (Miles et al. 2007), plus both 
positive and negative effects on skeleton formation 
(related to a decrease in carbonate saturation; Doney 
et al. 2009; Ries et al. 2009). 

The vulnerability of marine species and ecosystems to 
temperature, in particular, is well established (for reviews; 
Hoegh-Guldberg and Bruno 2010; Richardson et al. 2012; 
Wernberg et al. 2012). Conversely, the resilience of marine 
organisms to ocean acidification still remains a reasonably 
challenged concept (Dupont et al. 2010a; Hendriks and 
Duarte 2010). Recent meta-analyses assessing the biologi- 
cal effects of ocean acidification (Dupont et al. 2010a; 
Hendriks et al. 2010; Kroeker et al. 2010) concurred that 
it is unlikely to act in a uniform manner as variation exists 
in marine organism responses and resilience. Hence, if any 
meaningful comparisons are to be made on the response 
of marine organisms, they need to be hypothesis driven 
based on a priori assigned groupings, such as taxonomic 
groups or life stages (Dupont et al. 2010a). 

Studies of the biological effects of elevated temperature 
and acidification on marine organisms in isolation have 
provided some insight into the potential tolerance of spe- 
cies to these changing conditions (Gattuso and Hansson 
2009). However, given that these stressors are unlikely to 
operate independently (Halpern et al. 2007), there is now 
a need to gain a more ecologically realistic understanding 
of how the combined effects of temperature and acidifi- 
cation will affect marine biota (Sala et al. 2000; Fabry 
et al. 2008). This is vital in order to inform future adap- 
tive management strategies. Other recent meta-analyses, 
across ecological systems, have also shown that multiple 
stressors can lead to nonadditive interactions with 
responses dependent on the type of stressor as well as the 
level of ecological organization investigated (e.g., popula- 
tion vs. community, autotroph vs. heterotroph) (Crain 
et al. 2008; Darling and Cote' 2008; Tylianakis et al. 
2008). Moreover, the mechanism through which the stres- 
sor acts upon the organism will affect the response. Multi- 
ple stressors acting through a similar pathway may have an 
additive effect (Crain et al. 2008). In contrast, any stress- 
induced tolerances could lead to antagonisms (Blanck 
2002), while those stressors that act on different, but 
dependent mechanisms may act synergistically (Kneitel 
and Chase 2004). These reviews did, however, contain few, 



if any, studies that investigated both warming and acidifi- 
cation. Therefore, the concurrent effect of temperature 
and ocean acidification via elevated C0 2 remains unclear, 
but is likely to lead to complex biological outcomes. 

Organisms vary widely in their individual responses to 
ocean warming and acidification as a result of differences 
in their physiological and ecological characteristics 
(Dupont et al. 2008; Fabry 2008). For example, many 
marine organisms possessing a calcium carbonate 
(CaC0 3 ) structure would be considered more susceptible 
to ocean acidification as this process will impair their 
capacity to produce calcified skeletons (Doney et al. 
2009). Conversely, some species, including some calcified 
species, will have the capacity to buffer against the delete- 
rious effects of acidification by utilizing acid-base com- 
pensation (Claiborne and Evans 1992; Larsen et al. 1997), 
active mobility and metabolism (Widdicombe and Spicer 
2008; Whiteley 2011) or energy reallocation (Wood et al. 
2008; McDonald et al. 2009). Warmer temperature, up to 
a limit, stimulates metabolism in ectotherms, resulting in 
faster growth and development (Byrne 2011). Moreover, 
it has been speculated that warming could even amelio- 
rate the negative impacts of acidification (McNeil et al. 
2004; Kleypas and Yates 2009). 

Species responses to ocean warming and acidification 
are also likely to vary among life-history stages (Byrne 
2011). Early life-history stages are considered most suscep- 
tible to changes in both temperature and ocean acidification 
(Byrne 2011). These stressors may, however, have positive 
and/or negative effects for the successful recruitment of 
juveniles to the adult population. Trophic level is also 
likely to determine how species respond due to differences 
in environmental sensitivity (Petchey et al. 2004; Raffaelli 
2004). Previous study has suggested that the effects of 
multiple stressors are likely to act antagonistically in auto- 
trophs, but synergistically in heterotrophs (Crain et al. 
2008). Furthermore, as higher trophic levels contain less 
"biological insurance" (sensu Yachi and Loreau 1999), that 
is, less taxonomic, physiological, and genetic diversity, they 
are predicted to be more susceptible to multiple environ- 
mental perturbations (Christensen et al. 2006) which 
could act upon them synergistically (Crain et al. 2008). 

Using a meta-analytical approach of the peer-reviewed 
literature, we assessed the impacts and interactions of 
ocean acidification and warming on marine biological 
responses. Given that variability in the strength and direc- 
tion of responses was expected, we classified data accord- 
ing to taxonomic groups, calcifiers and noncalcifiers, life- 
history stage and level of trophic organization (autotroph 
and heterotroph) in terms of changes in rates of calcifica- 
tion, growth, photosynthesis, reproduction, and survival. 
Specifically, we aimed to address three questions: (i) How 
do marine organisms respond to warming and acidification 
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in isolation? (ii) How do marine organisms respond to 
the combined effects of warming and acidification? (iii) 
How do warming and acidification impacts interact? 

Material and Methods 

Data selection and suitability criteria 

Searches for peer-reviewed articles in which studies explic- 
itly investigated climate change using either elevated tem- 
perature, ocean acidification or elevated temperature and 
acidification were carried out using ISI Web of Science [v. 
5.8] © and Google Scholar using the following keywords: 
ocean warming, global warming, ocean acidification, 
hypercapnia, climate change, and combinations therein. In 
addition, we used the European Project on Ocean Acidifi- 
cation (EPOCA) blog (http://oceanacidification.wordpress. 
com/), citation searches; analysis of reference lists in 
comprehensive reviews (Hendriks et al. 2010; Kroeker 
et al. 2010; Wernberg et al. 2012), and then cross- 
referenced with the bibliographies of identified articles. 

We limited our analysis to studies published between 
1st January 1990 and 1st January 2012, as the majority of 
experimental climate change studies that manipulated cli- 
mate change conditions in line with IPCC AR1 predic- 
tions and subsequent updates (IPCC 1990, 2007) were 
published post 1990. Only controlled manipulative experi- 
ments were used for analysis. In addition, the control 
treatments of the environmental stressors (e.g., pH, C0 2 , 
or temperature) needed to represent current ambient lev- 
els and were based on the authors' opinion of "ambient". 
The experimental organisms had to be subjected to ele- 
vated temperature alone, acidification alone, or both 
warming and acidification. When studies included envi- 
ronmental variables in addition to temperature and ocean 
acidification (such as light availability or nutrients), these 
responses were only considered at "ambient" levels as 
determined by the authors'. To explore predicted future 
conditions for 2100, the manipulation treatments needed 
to conform to the IPCC IS92a "business-as-usual" emis- 
sion scenario for the year 2100 (IPCC 2007). We omitted 
studies that manipulated carbonate chemistry using acid 
addition, because it does not reproduce the changes in 
HC0 3 " concentration that occur as a result of increased 
C0 2 (aq) (Iglesias-Rodriguez et al. 2008a,b). Finally, only 
studies that reported a measureable biological response 
were included. 

As response variables, we used calcification (or dissolu- 
tion) rates, growth, photosynthesis, reproduction, and 
survival (mortality was converted to survival by using 1 - 
mortality). There were insufficient data on other response 
variables (e.g., feeding rates, metabolism) to allow quanti- 
tative analysis. A number of articles included more than 



one species, response, location, or treatment level and all 
were included in the analysis if they met the suitability 
criteria. This ensured that a broad range of responses 
could be fully explored, despite lessening the indepen- 
dence of the data from that particular study (Gurevitch 
et al. 1992). To maintain independence of data, we 
included only one response, chosen at random, from 
studies reporting several responses that could be classified 
in the same category (e.g., growth expressed as changes in 
length and biomass). Derived metrics from studies that 
included time-series data were based on the final time 
point of exposure. To investigate inherent biological vari- 
ability, records were categorized according to taxonomy, 
life-history stage, level of trophic organization (autotroph, 
heterotroph), and whether the organism possessed a 
CaC0 3 skeletal structure. 

To enable a calculation of effect size, studies that met 
our initial criteria could only be used if they reported a 
mean response value, some form of variance (standard 
deviation, standard error, or confidence interval), and a 
sample size. In some instances, data were extracted from 
graphical images in publications and in these situations 
data were extracted using the program GraphClick (v. 
3.0) (Neuchatel, Switzerland). 

Data analysis 

Biological responses to ocean warming and acidification 
were measured for each experiment to establish the pro- 
portional change between the control and treatment 
means using response ratios. In their original metric 
response ratios were weighted toward positive responses, 
so the response ratios were log transformed to maintain 
symmetry in the analysis and ease the biological interpre- 
tation (Hedges et al. 1999). We chose a log response ratio 
(InRR), over other methods, to estimate the effect size 
because of the high capacity to detect true effects and 
their robustness to small sample sizes (Lajeunesse and 
Forbes 2003). 

Analyses were carried out using the R (version 2.15.1; 
R Development Core Team 2012) package "Metafor" 
(Viechtbauer 2010). We selected a weighted random- 
effects model to estimate a summary effect size. Random- 
effects analysis assumes that the true effect size differs 
between experiments and the estimated summary effect is 
the mean of the effects observed across the studies. This 
meant that even if studies had a low weighting, the indi- 
vidual effect sizes from all of the studies could be incor- 
porated into the summary effect (Borenstein et al. 2009). 
This ensured that the biological variation inherent in the 
responses was properly accounted for. Both the within- 
study variance (inverse of the effect size variance) and the 
between-study variance (cr 2 p0 oied) were used to weight the 
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studies. Therefore, studies with higher replication and/or 
lower variance were considered more precise and 
weighted accordingly (Hedges and Olkin 1985). Between- 
study variance was estimated using the DerSimonian 
Laird method (DerSimonian and Laird 1986). 

Statistical significance was attributed to each summary 
effect size by calculating a bias-corrected 95% confidence 
interval (CI) (see Hedges and Olkin 1985) and comparing 
it with zero. If the summary effect size did not overlap 
zero then it was considered to be significantly different. A 
total heterogeneity statistic (Q) was used to ascertain that 
the variation observed was a combination of both true 
variation (between studies) and random error (within 
studies) (Borenstein et al. 2009). This was tested as the 
observed weighted sum of squares against a chi square 
distribution with n— 1 degrees of freedom, using the null 
hypothesis that observations shared a common effect size. 

Combinations of the treatment effect (C0 2 /pH, tem- 
perature, temperature and C0 2 /pH) and response vari- 
ables (calcification, growth, photosynthesis, reproduction, 
and survival) were used as the comparison groups in all 
analyses. Separate exploratory analyses were also used to 
test the differences between a priori defined groups using 
a mixed-effects model (see Viechtbauer 2010). It was 
appreciated that this form of multiple exploratory analy- 
ses on the same dataset could be prone to Type I error. 
Hence, we used these analyses to identify the underlying 
patterns of the biological responses. The categorical mod- 
erators used were the different taxonomic groups (corals, 
crustaceans, crustose coralline algae, echinoderms, fishes, 
noncalcifying algae, molluscs, phytoplankton, and seag- 
rasses), calcifying and noncalcifying organisms, develop- 
mental stages (embryos, larvae, juveniles, and adults), and 
trophic organization (autotroph and heterotroph). This 
process applied a summary effect size and 95% CI to each 
of the different categories for comparison. To formally 
test for differences between these categories, a test for het- 
erogeneity (Q M ) was used; this identified total heteroge- 
neity explained by that particular categorical moderator 
(Gurevitch et al. 1992). A significant Q M indicated that 
there was a difference between the categories. 

The taxonomic group of phytoplankton was initially 
divided into coccolithophores, cyanobacteria, diatoms, di- 
noflagellates, and foraminifera, however, results were 
pooled after detecting no differences using a test for het- 
erogeneity (Q M )- When assigning a "trophic organization" 
to each observation, we defined autotrophs as any organ- 
ism capable of producing organic carbon-based com- 
pounds from inorganic sources, through either 
photosynthesis or chemosynthesis. For example, corals 
were designated as autotrophic as photosynthesis is gener- 
ally regarded as their principal mode of carbon acquisi- 
tion (Hoogenboom et al. 2006; Mass et al. 2007). Over all 



of the meta-analytical results, the summary effect sizes 
were not reported if there were fewer than five studies 
available for analysis, and categorical moderators were 
not reported if there were fewer than three studies. This 
was a pragmatic decision to ensure that a broad range of 
responses could be assessed, as some categories only had 
a few studies that met our criteria. Therefore, the categor- 
ical analyses did not always include all the observations 
from the full model. 

Interactions between multiple stressors 

Interaction strength between ocean warming and acidifi- 
cation was ascertained according to the methods for fac- 
torial meta-analysis (Gurevitch et al. 2000; Crain et al. 
2008). To be included, studies needed to be controlled 
factorial experiments reporting four outcomes of acidifi- 
cation [Yacidli warming [F temp ], acidification and warming 
[it>oth]> an d a control treatment [Y^] (Underwood 1997). 
Therefore, not all the observations from the full model 
could be analyzed. Multiple observations from the same 
study were included if separate factorial results were pro- 
vided. The interaction strength (lnRR int ) and individual 
effects (lnRR ac id and lnRR, emp ) of each study were then 
calculated as 



lnRR ir 



[ temp ) 



2s 



InRR 



acid 



(^acid 



InRR, 



temp 



where Y is the biological response for the indicated group 
(in subscript), and s is the pooled standard deviation. The 
sampling variance of lnRR int is 



s^(lnRR int ) 



-Nadd + Ntemp + Nboth + N< 



1 



+ ■ 



1 



- + - 



1 



+ 



lnRR int 

!(N acid + N temp + Nboth + ' 



where N is the sample size for the indicated group (in 
subscript) . 

Although less conservative than an additive model (Folt 
et al. 1999), we used a multiplicative model to calculate 
the interactions as the underlying model of the metric 
InRR is multiplicative (Hawkes and Sullivan 2001; Morris 
et al. 2007), and this model is also thought to be more 
biologically realistic (Sih et al. 1998). Synergisms occur 
when the cumulative response of a stressor 
(iboth — itemp) is greater than the stressor in isolation 
(^add — Y ct ), and antagonistic when the cumulative 
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impact is smaller than expected (Folt et al. 1999). The 
interaction strength (lnRR int ) was then synthesized into a 
summary effect for each type of biological responses using 
a random-effects model with the same methodology as 
outlined previously. Therefore, the interaction type was 
classified as multiplicative (i.e., the null hypothesis) if the 
interaction effect size and 95% CI overlapped zero. If the 
direction of both individual stressors was positive, the 
interaction was considered synergistic when the interac- 
tion effect size was greater than zero, and antagonistic if 
less than zero. Similarly, when the direction of the indi- 
vidual stressors was either both negative or had opposite 
signs, then the interpretation was reversed (i.e., synergies 
occurred when the interaction effect size was less than 
zero). 

Sensitivity analyses and publication bias 

Sensitivity analysis was used to investigate the influence 
of any experimental study that demonstrated an unusually 
large effect size. This was achieved in a stepwise manner 
by ranking each experiment by the magnitude of effect 
size, removing the largest one, and rerunning the analysis. 
Likewise, if any study contributed five or more observa- 
tions to a category, the study was omitted and the analy- 
ses were rerun. If studies were considered to be driving 
the results, then they were omitted from the analysis of 
that response variable. 

The number of studies with an effect size of zero that 
would be required to change the results of the meta anal- 
ysis from significant to nonsignificant ("file drawer prob- 
lem") was determined using Rosenthal's fail-safe number 
(Rosenthal 1979). It was decided that if five or less studies 
(of zero effect size) were required to change the effect 
size, then that categorical analysis was not considered 
robust. 

Results 

Overall biological responses 

Out of 196 peer-reviewed articles that investigated the 
biological responses of marine organisms to ocean warm- 
ing and/or acidification 107 met our criteria, giving 623 
unique observations (Table SI). Observations that did not 
meet the selection criteria are listed in Table S2, and the 
results from all the heterogeneity tests for overall within 
effects (Q) and between categories (Q M ) are reported in 
Table S3. 

Meta-analysis of the whole dataset revealed that calcifi- 
cation and reproduction were negatively affected by ocean 
acidification and neutrally affected by ocean warming. In 
contrast, the independent effects of ocean acidification 



and warming had no effect on growth or photosynthesis, 
while both ocean acidification and warming had signifi- 
cant negative effects on survival (Fig. 1). 

Combined warming and acidification had a significant 
negative effect on calcification, survival, and reproduction 
with the magnitude of response greater than that 
observed for the stressors in isolation (Fig. 1). Concurrent 
acidification and warming had no effect on growth over- 
all, but had a significant positive effect on photosynthesis 
in autotrophs (Fig. 1). 

Significant within-group heterogeneity was observed in 
all the responses (Table S3), therefore in order to quantify 
patterns in the literature regarding the effects of warming 
and acidification, we made comparisons between the a 
priori groupings (taxonomic groups, calcifiers and non- 
calcifiers, life-history stage, and level of trophic organiza- 
tion) for each biological response. 

Calcification 

Although the independent effects of ocean acidification 
and warming on calcification were varied, we did not 
detect differences between taxonomic groups (Fig. 1; 
Table S3), life-history stages (Fig. 3; Table S3), or trophic 
organization (Fig. 4; Table S3). Conversely, the combined 
effects of ocean warming and acidification on calcification 
varied significantly by taxonomic group (Fig. 1; 
P < 0.001; Table S3), with echinoderms more negatively 
affected than either corals or crustaceans (both 
P < 0.001); life-history stage (Fig. 3; P = 0.0397; Table 
S3) with juveniles more negatively affected compared with 
adults (P = 0.040); and trophic organization (Fig. 4; 
P = 0.010; Table S3) with heterotrophs being negatively 
affected by warming and acidification, while autotrophs 
were neutrally affected (Fig. 4). 

Growth 

Despite overall growth responses being neutrally affected 
by ocean acidification and warming, we detected differ- 
ences in the combined effects between taxonomic groups 
(Fig. 1; P < 0.001; Table S3). Crustaceans and molluscs 
were both significantly negatively affected, while echino- 
derms and phytoplankton were positively affected 
(Fig. 1). The combined effects of warming and acidificat- 
ion varied with trophic organization (Fig. 4; P < 0.001; 
Table S3), with autotrophs being positively affected and 
heterotrophs negatively affected (Fig. 4). Similarly, ocean 
warming also varied between taxonomic groups (Fig. 1; 
P < 0.001; Table S3), with echinoderms, phytoplankton, 
fishes, and macroalgae positively affected, and molluscs 
negatively affected (Fig. 1); and trophic organization 
(Fig. 4; P = 0.004; Table S3) with a positive effect on 
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Figure 1. The mean effect of ocean acidification (clear circles), ocean warming (gray circles), and combined ocean acidification and warming 
(black circles) on calcification, growth, photosynthesis, reproduction, and survival for different taxonomic groups. The mean log response ratio 
and ±95% confidence intervals are shown for overall (combined results), calcifiers (calcifying algae, corals, crustaceans, echinoderms, molluscs, 
and phytoplankton) and noncalcifiers (fishes, noncalcified algae, seagrass). The number of observations in each analysis is shown in parentheses. 
The zero line indicates no effect, and significance (*) of mean effects is determined when the ±95% confidence interval does not overlap zero. 



autotrophs and a neutral effect on heterotrophs. Ocean 
warming also positively affected noncalcifiers (Fig. 2). 
The independent effects of ocean acidification varied by 



taxonomic group (Fig. 1; P < 0.001; Table S3), with crus- 
tose coralline algae, corals, molluscs, and phytoplankton 
negatively affected, crustaceans and echinoderms unaffected, 
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Figure 2. The mean effect of ocean acidification (clear circles), ocean warming (gray circles), and combined ocean acidification and warming 
(black circles) on growth, photosynthesis, and survival for calcifying and noncalcifying organisms. The mean log response ratio and ±95% 
confidence intervals are shown for calcifiers and noncalcifiers. The number of observations in each analysis is shown in parentheses. The zero line 
indicates no effect, and significance (*) of mean effects is determined when the ±95% confidence interval does not overlap zero. 



and fishes, macroalgae, and seagrass positively 
affected (fishes and macroalgae P < 0.05, seagrass 
P = 0.06); and calcifiers/noncalcifiers (Fig. 2; P < 0.001; 
Table S3) with the calcifying organisms exhibiting gener- 
ally negative effects and noncalcifiers generally positive 
effects (Figs. 1, 2). No differences in growth responses 
were detected between life stages with any of the stressors 
(Fig. 3). 

Photosynthesis 

We did not detect differences in the effects of combined 
warming and acidification on photosynthesis between 
either taxonomic groups (corals and phytoplankton; 
Fig. 1), calcifiers/noncalcifiers (Fig. 2), or life stages 
(Fig. 3), despite finding significant differences between 
these a priori groupings in both acidification and warm- 
ing independently (Table S3). We did, however, find that 
in combination significant positive effects were observed 
in phytoplankton (Fig. 1) and calcifying organisms 
(Fig. 2). Similarly, ocean warming positively affected pho- 
tosynthesis in phytoplankton and calcifying organisms, 
but had a negative effect on both macroalgae (Fig. 1) and 
noncalcifying organisms (Fig. 2). The independent effects 
of ocean acidification on photosynthesis also varied by 
taxonomic group (Fig. 1; P < 0.001; Table S3), with cor- 
als being negatively affected, and both phytoplankton and 
seagrass positively affected. Ocean acidification also posi- 



tively affected photosynthesis in noncalcifying organisms 
(Fig. 2). 

Reproduction 

The effects of warming and acidification, both combined 
and independently, on reproduction were varied between 
the taxonomic groups (Fig. 1; all P < 0.001; Table S3). In 
response to ocean warming alone positive responses by 
crustaceans and macroalgae and autotrophs were observed 
while heterotrophs responded negatively (P = 0.0071, 
Table S3). In all three combinations of stressors, echino- 
derms were unaffected while molluscs were negatively 
affected. 

Survival 

The combined effects of warming and acidification 
resulted in significant negative responses across the cate- 
gories, but only varied by taxonomic group (Fig. 1; 
P = 0.030; Table S3). Specifically, combined acidification 
and warming had a larger negative effect on molluscs 
than crustaceans (Fig. 1; P = 0.029). Independently, the 
effects of ocean warming on survival also varied by life 
stage (Fig. 3; P < 0.001; Table S3) with a larger negative 
effect on juveniles compared with adults (P = 0.023) and 
on larvae compared with both juveniles and adults (both 
P < 0.001). Additionally, ocean warming also had a 
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Figure 3. The mean effect of ocean acidification (clear circles), ocean warming (gray circles), and combined ocean acidification and warming 
(black circles) on calcification, growth, and survival in different life-history stages. The mean log response ratio and ±95% confidence intervals are 
shown for embryos, larvae, juveniles, and adults. The number of observations in each analysis is shown in parentheses. The zero line indicates no 
effect, and significance (*) of mean effects is determined when the ±95% confidence interval does not overlap zero. 



negative effect on molluscs (Fig. 1), calcifiers (Fig. 2), and 
both autotrophs and heterotrophs (Fig. 4). Survival 
responses to ocean acidification in isolation did not 
vary between taxonomic group, calcifiers/noncalcifiers, 
life stage, or trophic organization. However, ocean 
acidification did negatively affect corals and molluscs 
(Fig. 1), calcifiers (Fig. 2), larvae (Fig. 3), and autotrophs 
(Fig. 4). 

Interactions between multiple stressors 

Interaction strength between ocean warming and acidifi- 
cation was determined for the subset of fully factorial 
studies, and therefore did not include all the observations 
from the full model. Significant synergistic interactions 
were observed for calcification (z = 3.69, P < 0.001), pho- 
tosynthesis (z = 2.23, P = 0.026), reproduction (z = 2.97, 
P = 0.003), and survival (z = 3.04, P = 0.002), but not 
for growth (z = 1.56, P = 0.117). 

Sensitivity analyses and publication bias 

To test the robustness of our analyses against large effect 
sizes, we removed each comparison stepwise and reran 
each analysis, omitting experiments if they changed the 
significance of either heterogeneity or the mean effect size 
of the response variables. This resulted in 12 experiments 
being omitted from subsequent analyses across several 
treatment-response variable scenarios (see Table S2 for 



more detail). We used Rosenthal's fail-safe number to 
assess the importance of potential publication bias and 
found that our response variables were robust, with the 
lowest value being 358 additional studies required to 
change the effect size (based on an original experiment 
quantity of 11). No individual study contributing more 
than five experiments changed the significance of either 
the heterogeneity or mean effect size of the response vari- 
ables. 

Discussion 

This study used a meta-analytical approach to assess the 
impacts and interactions of ocean acidification and warm- 
ing on marine biological responses. Responses were classi- 
fied according to taxonomic groups, calcifiers and 
noncalcifiers, level of trophic organization (autotroph and 
heterotroph), and life-history stage in terms of changes in 
rates of calcification, growth, photosynthesis, reproduc- 
tion, and survival. We assessed marine organisms 
responses to warming and acidification in isolation and 
in combination and assessed the biological implications 
of their interaction. 

Effects of ocean acidification 

Ocean acidification generally had an adverse effect on a 
large range of marine biota with more specific differ- 
ences between life-history characteristics. Calcifying 
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organisms were generally negatively affected by ocean 
acidification, with significant biological variability in 
responses, such as molluscs being negatively affected, 
which may be due to their poor ion regulation and 
inability to buffer their internal compartments (Fabry 
et al. 2008; Widdicombe and Spicer 2008; Melzner et al. 
2009; Dupont et al. 2010a). Conversely, crustaceans 
were generally unaffected by acidification perhaps 
due to their active mobility, higher metabolism, 
and capacity to control intracellular pH (Gaillard and 
Malan 1983; Widdicombe and Spicer 2008; Whiteley 
2011). 

Overall, less is known about how noncalcifying organ- 
isms are likely to respond to acidification (Connell and 
Russell 2010), particularly marine fishes (Ishmatsu et al. 
2008; Munday et al. 2009), however, our results show 
that noncalcifiers were generally unaffected by acidificat- 
ion and for growth were positively affected when analyzed 
together as a group. Fishes are thought to have more effi- 
cient acid-base regulation compared to invertebrates 
(Widdicombe and Spicer 2008), which, when coupled 
with an increased food intake or reduced energy expendi- 
ture (Munday et al. 2009), could explain the positive 
growth response observed. Similarly, noncalcifying marine 
autotrophs demonstrated an increased growth to ocean 
acidification. This is possibly via their capacity to derive 
dissolved inorganic carbon from the increased C0 2 (aq) 
(Beer and Koch 1996). 



Autotrophs may be capable of increasing their inor- 
ganic carbon assimilation (Rost and Riebesell 2004) and 
buffering the negative effects on calcification (Ries et al. 

2009) . However, negative effects on growth were observed 
in calcifying autotrophs. Our results suggest that the neg- 
ative effects on calcifying autotrophs may outweigh posi- 
tive effects associated with the increased availability of 
C0 2 (aq) as a substrate for photosynthesis. We also show 
some evidence that calcification and survival in early life- 
history stages were more negatively affected by ocean 
acidification, highlighting not only the susceptibility of 
early life stages, but also the subtle nature of life-history 
responses when compared to species-specific effects of 
heterogeneity (Kurihara 2008; Kroeker et al. 2010). 
Importantly, our results agree with the findings from pre- 
vious meta-analyses (Dupont et al. 2010a; Kroeker et al. 

2010) while introducing a further 48 key studies (168 
additional data points). 

Effects of ocean warming 

Moderate elevations in temperatures will increase meta- 
bolic rates (Hochachka and Somero 2002), which influ- 
ences key biological processes that regulate life-history 
characteristics (O'Connor et al. 2007). While marine 
organisms are capable of acclimation to a range of tem- 
peratures, once their thermotolerance limits are exceeded, 
organism fitness is reduced and the risk of mortality 
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increases (Hofmann and Todgham 2010; Tomanek 2010). 
We found that both noncalcifying organisms and auto- 
trophs demonstrated increased growth under warming 
conditions, likely due to this increase in metabolic rate 
(Hochachka and Somero 2002). Warming, however, had 
no effect on the growth of heterotrophs and even nega- 
tively affected their calcification. Possibly because the 
metabolism complexes of autotrophs (photosynthesis-lim- 
ited) are less sensitive to ocean warming than the respira- 
tion-limited metabolism of heterotrophs (Lopez-Urrutia 
et al. 2006). Moreover, our results also support the 
hypothesis that the threshold for deleterious warming 
may vary between developmental stages (Byrne et al. 
2009, 2010) with survival significantly lower in larvae 
than in juveniles, and significantly lower in juveniles than 
in adults. 

Simultaneous acidification and warming 

Meta-analysis of the full dataset revealed that the com- 
bined stressors caused significant negative effects on calci- 
fication, reproduction, and survival, and a significant 
positive effect on photosynthesis, but no effect on growth. 
Importantly, we also found that four of the five responses 
(calcification, photosynthesis, reproduction, and survival) 
showed a synergistic interaction between acidification and 
warming. Although, such synergistic interactions between 
stressors (i.e., where the outcome was greater than the 
sum of the individual stressors; Folt et al. 1999) are rela- 
tively common (e.g., Sala et al. 2000; Harley et al. 2006), 
they are concerning because they are also unpredictable. 
Hence, such synergies limit our capacity to predict poten- 
tial future impacts from single-stressors studies. 

Ecological synergies are important to marine systems 
(Paine et al. 1998; Harley et al. 2006; Sutherland et al. 

2006) because they can further exacerbate adverse effects 
and reduce ecosystem resilience (Folke et al. 2004). They 
can also introduce indirect effects via biotic interactions 
(Darling and Cote' 2008; Tylianakis et al. 2008). For 
example, climate-driven changes in plankton communities 
can regulate top predators through bottom-up control 
(Beaugrand et al. 2003). Moreover, as marine systems are 
subject to multiple interacting stressors (Halpern et al. 

2007) , it is possible that the addition of further stressors 
would introduce additional adverse consequences (e.g., 
Przeslawski et al. 2005). Therefore, our results highlight 
the need to move away from single-stressor studies and 
toward more ecologically realistic research incorporating 
multiple stressors, in order to more fully understand how 
near-future anthropogenic change will affect marine bio- 
diversity. 

There was, as expected, variation in our analyses 
among the biological responses to combined warming 



and acidification between different taxonomic groups, cal- 
cifiers and noncalcifiers, trophic levels, and life-history 
stages. The combination of warming and acidification 
generally exhibited a stronger effect (either positive or 
negative) than when exposed to the stressors in isolation. 
For instance, echinoderms are highly vulnerable to ocean 
acidification (Dupont et al. 2010b), likely due to their 
skeletons being formed from highly soluble magnesium 
calcite (Politi et al. 2004). However, the addition of mod- 
erate increases in temperature (+4°C), as predicted for 
the end of the century, resulted in further adverse effects 
(e.g., Byrne et al. 2011), such as the highly negative calci- 
fication responses shown here. In some instances, the 
combined effects of warming and acidification were 
reduced compared with the individual effects. Corals had 
both calcification and photosynthesis negatively affected 
by ocean acidification in isolation, while they were unaf- 
fected by the combined effects, suggesting that the addi- 
tion of warming may ameliorate the adverse effects of 
acidification (McNeil et al. 2004; Kleypas and Yates 
2009). These differences in the resilience of marine organ- 
isms will have important implications for ecosystem level 
responses. 

Interestingly, the combined effect of warming and acid- 
ification positively affected growth in echinoderms, which 
could be explained by energy allocation, where the cost of 
homeostatic regulation can be influenced by changes in 
somatic and reproductive growth performance (Melzner 
et al. 2009). Hence, if more energy is utilized to maintain 
growth, then calcification responses could be more 
adversely affected (e.g., Arnold et al. 2009; McDonald 
et al. 2009). Other studies have shown an alternative 
strategy for energy allocation where growth was negatively 
affected, while calcification was maintained (e.g., Wood 
et al. 2008; or crustaceans, this study), highlighting the 
species-specific nature of biological responses. 

Heterotroph responses to ocean acidification and 
warming individually differed, with acidification reducing 
growth, but not affecting calcification or survival; while 
warming alternatively did not affect growth, reducing 
both calcification and survival responses. However, the 
energetic demands of the combined effects of acidification 
and warming on heterotrophs resulted in calcification, 
growth, and survival all being reduced. Conversely, the 
combined effects of warming and acidification positively 
affected growth in autotrophs, likely due to the effect of 
temperature on metabolic rate (Hochachka and Somero 
2002), while C0 2 , acted as a substrate for photosynthesis 
and possibly indirectly promoted growth (e.g., phyto- 
plankton; Loehle 1995). Moreover, as multiple stressors 
affecting autotrophs are likely to act antagonistically 
(Crain et al. 2008), and the dissolved inorganic carbon 
sources utilized by marine autotrophs is set to increase 
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(Raven 2005), photosynthesizing organisms will likely be 
more resilient to conditions predicted for the end of the 
century; as long as they do not exceed their thermotoler- 
ance (Hofmann and Todgham 2010; Tomanek 2010) and 
are not limited by other factors, such as inorganic nutri- 
ent availability (Langdon and Atkinson 2005; Cohen and 
Holcomb 2009; Ries et al. 2009). 

Our results also show that the combined effects of 
warming and acidification had significant negative effects 
on the survival of early life-history stages. Although there 
was insufficient data to compare across all life-history 
stages, both larvae and juveniles were highly susceptible 
to changes in temperature and ocean acidification, sup- 
ported by previous research (Gosselin and Qian 1997; 
Hunt and Scheibling 1997; Byrne 2011). Indirectly, ocean 
warming can reduce the mortality of larvae by shortening 
the planktonic duration (Lamare and Barker 1999), when 
they are most vulnerable to predation (O'Connor et al. 
2007). Our results, however, identified that the combined 
effects of ocean acidification and warming increased mor- 
tality, indicating that multiple stressors will have impor- 
tant implications for population persistence, potentially 
acting as a bottleneck for some species (Dupont et al. 
2010b; Byrne 2011). 

Conclusions 

Quantitative syntheses of the published literature can pro- 
vide powerful inferences, however, like all analyses they 
are subject to caveats. We identified and incorporated the 
available literature that met our selection criteria, how- 
ever, this also outlines the current gaps in knowledge and 
highlights opportunities for further study. Moreover, spe- 
cies-specific sources of heterogeneity are always likely to 
make some results from the literature greatly context 
dependent (e.g., Fabry 2008; Kurihara 2008; Dupont et al. 
2010a; Hendriks et al. 2010; Kroeker et al. 2010). Hence, 
our findings highlight the complexity of marine organism 
responses to ocean warming, acidification, and their inter- 
action. The magnitude, direction, and interaction of the 
effects varies between response types, likely a result of the 
pathways driving the biological responses. Responses also 
differ between taxonomic groups, trophic levels, and life- 
history stages. Most importantly, we observed synergistic 
interactions between ocean acidification and warming in 
four of the five biological responses measured (calcifica- 
tion, photosynthesis, reproduction, and survival), 
highlighting the difficulties in making inferences from 
single-stressor studies. However, single-factor studies in 
junction with those that manipulate multiple stressors can 
play a vital role in understanding the pathways through 
which particular stressors operate and will enable a more 
accurate assessment of the likely outcomes of interactions 



between warming and acidification. Importantly, we must 
also consider further abiotic and biotic stressors in the 
marine environment that are likely to also interact with 
warming and acidification (Halpern et al. 2007) as well as 
scaling up studies from individuals and populations to 
communities and ecosystems (Harley et al. 2006). Such 
large-scale multifactorial experiments would not only 
increase our knowledge of the functioning and resilience 
of marine ecosystems, but provide explicit evidence to 
policymakers on the effectiveness of conservation and 
management strategies in response to future environmen- 
tal change. 
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